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Abstract. The solar wind conditions at one astronomical unit (AU) can be strongly 
disturbed by the interplanetary coronal mass ejections (ICMEs). A subset, called magnetic 
clouds (MCs), is formed by twisted flux ropes that transport an important amount of 
magnetic flux and helicity which is released in CMEs. At 1 AU from the Sun, the magnetic 
structure of MCs is generally modeled neglecting their expansion during the spacecraft 
crossing. However, in some cases, MCs present a signiflcant expansion. We present here 
an analysis of the huge and significantly expanding MC observed by the Wind spacecraft 
during 9 and 10 November, 2004. This MC was embedded in an ICME. After determining 
an approximated orientation for the flux rope using the minimum variance method, we 
precise the orientation of the cloud axis relating its front and rear magnetic discontinuities 
using a direct method. This method takes into account the conservation of the azimuthal 
magnetic flux between the in- and out-bound branches, and is valid for a flnite impact 
parameter (i.e., not necessarily a small distance between the spacecraft trajectory and 
the cloud axis). The MC is also studied using dynamic models with isotropic expansion. 
We have found (6.2 ± 1.5) x 10^° Mx for the axial flux, and (78 ± 18) x 10^° Mx for the 
azimuthal flux. Moreover, using the direct method, we find that the ICME is formed by 
a fiux rope (MC) followed by an extended coherent magnetic region. These observations 
are interpreted considering the existence of a previous larger flux rope, which partially 
reconnected with its environment in the front. We estimate that the reconnection process 
started close to the Sun. These findings imply that the ejected flux rope is progressively 
peeled by reconnection and transformed to the observed ICME (with a remnant flux rope 
in the front part). 

Keywords: Coronal Mass Ejections, Interplanetary; Magnetic fields. Interplanetary; Mag- 
netic Reconnection, Observational Signatures; Solar Wind, Disturbances 



1. Introduction 



A magnetic configuration, previously in equilibrium in the solar atmosphere, 
can reach a global instability threshold when the magnetic stress becomes 
too high. In this case, the plasma is ejected into the interplanetary (IP) 
medium and is observed as a Coronal Mass Ejection (CME) by the solar 
coronagraphs. This magnetized mass, which can be expelled as fast as few 
times 1000 km/s, is recognized in the IP space as an interplanetary CME, 
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ICME, see e.g., Wimmer-Schweingruber et al. (2006). During its travel from 

the Sun to 1 AU, fast CMEs are slowed down due to drag forces between the 
ICME and the solar wind environment (see, e.g., Vrsnak and Gopalswamy, 

2002) . Thus, at 1 AU they can reach speeds as high as ~ 1000 km/s. 

1.1. Magnetic Clouds 

Magnetic clouds (MCs) are a particular subset of ICMEs. They are formed 
by twisted magnetic flux tubes that carry a large amount of magnetic helicity 
from the Sun to the IP medium. They also transport significant amounts of 
magnetic flux, mass, and energy. The principal characteristics of these mag- 
netic structures are: (i) an enhanced magnetic field, (ii) a smooth rotation 
of the magnetic field vector through a large angle (r^ 180°), and (iii) a low 
proton temperature (Klein and Burlaga, 1982). 

The magnetic field in MCs is relatively well modeled by the so-called 
Lundquist's model (Lundquist, 1950), which considers a static and axially- 
symmetric linear force-free magnetic configuration (e.g., Goldstein, 1983; 
Burlaga, 1988; Lepping et al, 1990; Burlaga, 1995; Lynch et al, 2003). 
However, many other different models have been also used to describe the 
magnetic structure of MCs. 

Some models consider the MC as a rigid body during the time it travels 
through the solar wind and crosses the spacecraft. Farrugia et al. (1999) 
considered a cylindrical shape for the cloud cross-section and a non-linear 
force-free field; while Mulligan et al. (1999), Hidalgo et al. (2002), and 
Cid et al. (2002) considered a cylindrical cloud but a non-force free field. 
Non cylindrical static models have been also applied to MCs (e.g., Hu and 
Sonnerup, 2001, Vandas and Romashets, 2002). A comparison of global 
quantities (magnetic fluxes and helicity) derived from different static models 
has been done by Dasso et al. (2003) and Dasso et al. (2005b). Different tech- 
niques have been compared using synthetic data and analyzing the output 
of numerical simulations of MCs (Riley et al., 2004). 

Some MCs present a signiflcantly larger velocity in their front than in 
their back, a characteristic of expansion. Thus, some authors have used 
dynamical models to describe these clouds during the observation time; 
they have considered two cases: only with a radial expansion (see, e.g., 
Farrugia et al, 1993, Osherovich et al, 1993, Farrugia et al., 1997, Nakwacki 
et al., 2005), and with expansion in the radial and axial directions (see, e.g., 
Shimazu and Vandas, 2002, Berdichevsky et al, 2003). Some dynamical 
models consider an expanding elliptical shape for the MC (e.g., Hidalgo, 

2003) . The main aim of these models is to take into account the evolution 
of the magnetic fleld as the spacecraft crosses the MC; then, to correct the 
effect of mixing spatial- variation/time-evolution in the observations to get 
a better determination of the MC field (and related characteristics). 
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1.2. Aims of this study 

We analyze the MC detected inside the ICME observed at LI between 
Nov. 9, 2004, at 20:25 UT and Nov. 11, 2004, at 18:45 UT (Harra et al, 
2006). This is a very fast, left-handed and huge MC with a size larger than 
0.2 AU in the Earth-Sun direction. It has a very intense magnetic field 
(> 40 uT), and expands strongly, with a velocity of ~ 850 km/s in its front 
and ~ 600 — 700 km/s in its back, depending on where the rear boundary is 
set. This MC presents one of the largest (ever observed) velocity differences 
between its front and its back (Nakwacki et al, 2007). 

We analyze this MC using a model-independent method, called direct 
method (Dasso et al., 2006). It takes into account the magnetic flux con- 
servation in closed structures, such as flux ropes. This method gives us an 
estimation of the magnetic flux in the MC directly from the data, and allows 
us to improve the determination of the orientation of the MC axis, as well 
as its boundaries. Finding the boundaries for some MCs is an open issue 
(Russell and Shinde, 2005; Wimmer-Schweingruber et al, 2006). Indeed, 
several authors, using different proxies, set them at different times with the 
consequent differences in the estimated MC axial orientation and in the 
estimations of global magnctohydrodynamic quantities. 

In Section 2, we describe the magnetic and plasma properties of the MC. 
Then, in Section 3, we present the direct method and the results obtained. 
In Section 4, wc fit the two dynamical models to velocity and magnetic field 
observations. From the fitted models wc compute the magnetic fluxes, and 
compare them with the values obtained from the static Lundquist's model 
and from the direct method. Finally, in Section 5, we give a summary and 
our conclusions. 



2. The studied cloud 

In this section we analyze in situ magnetic and plasma observations of the 
MC observed on Nov. 9-10, 2004, which is located inside an ICME. We define 
a local frame of coordinates, attached to the MC, and we analyze the data 
in this frame. 

2.1. One or two clouds? 

An ICME was observed from 20:25 UT on 9 Nov. to 18:45 UT on 11 Nov. 
The ICME is preceded by a typical piled up solar wind material (panel Hp 
in Figure 1). This corresponds to plasma and magnetic field pushed from 
behind by the ICME, forming the turbulent prc-ICME sheath (notice the 
high level of fluctuations in 9b)- A forward shock is located in front of the 
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sheath (at 9:15 UT of Nov. 9). This event presents also a non-typical second 
shock inside the sheath that precedes the ICME (at 18:20 UT of Nov. 9). 
For a deeper description of shocks and their association with solar sources 
see Harra et al. (2006), in particular their Table I. 

The end of the ICME is marked by a thick solid line (Nov. 11, 18:45 UT, 
label "end" in Figure 1) as defined by Harra et al. (2006). After Nov. 11 at 
18:45 UT, the magnetic field is consistent with Parker's spiral (changes of 
magnetic sectors are observed in 0^). The higher level of magnetic fluctua- 
tions, typical of the fast solar wind, confirms this interpretation. 

Two MCs were initially reported inside this ICME: the first one starting 
on Nov. 09, at 20:54 UT, and finishing on Nov. 10, at 03:24 UT, and the 
second one starting on Nov. 10, at 03:36 UT, and finishing on Nov. 10 
at 11:06 UT (http://lepmfi.gsfc.nasa.gov/mfi/mag_cloud_Sl- 
.html). The rotation of the field is indeed larger than usual (close to one 
turn. Figure 1). Moreover, the magnetic field observed has indeed a non- 
classical structure: it is very strong in the front, progressively decreases in 
the back and has an extended weak tail. Still, a coherent and continuous 
variation of both the field strength and the plasma velocity is observed, 
without clear evidence of two independent magnetic structures. 

Other studies concluded that only one extended MC was present. In Harra 
et al. (2006) and Longcope et al. (2007) the boundaries of the magnetic cloud 
were selected as starting on Nov. 9, at 20:30 UT, and finishing on Nov. 10, 
at 10:00 UT. Similar boundaries (starting on Nov. 9, at 20:40 UT, and 
finishing on Nov. 10, at 10:20 UT) were chosen by Qiu et al. (2007). Both 
the direct method and data modeling confirm the presence of only one fiux 
rope (Sections 3 and 4). 

2.2. Summary of the cloud observations 

We analyze the in situ measurements of the magnetic field components 
obtained by the Magnetic Field Instrument (MFI, Lepping et al, 1995) 
and the plasma quantities obtained by the Solar Wind Experiment (SWE, 
Ogilvie et al, 1995), both aboard Wind. There is a small data gap in SWE 
from Nov. 9, 21:58 UT to 22:27 UT (see Figure 1). 

The magnetic field observations are in GSE (Geocentric Solar Ecliptic) 
coordinates. In this right-handed system of coordinates, xgse corresponds 
to the Earth-Sun direction, zqse points to the North (perpendicular to the 
ecliptic plane) and yosE is in the ecliptic plane and points to the dusk when 
an observer is near Earth (thus, opposing to the planetary motion). 

The front boundary of the MC (Nov. 9, at 20:30 UT) is well defined 
(vertical dashed line in Figure 1, label "in"). The magnetic field presents a 
North- West-South rotation with time (see 9b and (pB panels in Figure 1); 
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Figure 1. Wind observations for the magnetic cloud observed inside the ICME of 9-11 
Nov. 2004 (time cadence of 100 seconds). Prom upper to lower panels: absolute value of 
the magnetic field {B = \B\), latitude {9b) and longitude {6b) angles of B in GSE, the 
bulk velocity (V^), the proton density (rip), and the proton plasma beta (/3p), all of them 
as a function of time. The vertical dashed line "in" corresponds to the start of the field 
coherent rotation (beginning of the MC, Nov. 9 at 20:30 UT), the vertical dotted line to 
the cloud center (Nov. 10 at 02:02 UT, see Section 3.3), the vertical dash-dotted lines 
"outi" and "outs" to the extremes of the range of possible endings for the rotation of B 
(end of the cloud, Nov. 10 at 08:15 UT and 10:20 UT), and the dash-dotted line "out2" 
to a strong discontinuity (in between this range, on Nov. 10, at 09:00 UT). The thin solid 
vertical line "back" marks a strong magnetic discontinuity (Nov. 11, at 06:02 UT) and the 
thick vertical sohd line "end" the end of the ICME (Nov. 11, at 18:45 UT). Horizontal 
dotted lines in 6b and <j)B panels indicate the orientation of Parker's spiral, while the one 
in (5p panel marks the mean value for a set of MCs studied by Lepping et al. (2003). 



thus, the MC is formed by a left-handed flux rope with its axis almost on 
the ecliptic and pointing to the west {yaSE < or ~ 270°). 

A characteristic of this MC is its very strong expansion (see panel V in 
Figure 1). The observed plasma velocity, V, goes from ~ 850 km/s at the 
beginning to ~ 700 km/s close to the cloud end, a difference of 150 km/s 
in the observed time range (~ 15 hours). This implies an expansion of ~ 
10 km/s per hour. 
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The MC rear boundary is uncertain, the vertical dash-doted hues in 
Figure 1 indicate a possible range (from Nov. 10 at 8:15 UT to 10:20 UT, 

labels "outi" and "outs", respectively). However, the decrease of \B\ and 
the strong expansion are still present at later times indicating a backward 
extension of the MC. The mean value of f3p (ratio of the proton pressure 
to the magnetic pressure) is 0.12 for a sample of MCs studied by Lepping 
et al. (2003) (horizontal dotted line in panel /3p of Figure 1). Then, another 
indication that the MC is more extended in the back is the presence of 
Pp < 0.12 after 10:20 UT. This region extends up to a strong discontinuity 
in Ob and density (on Nov. 11 at 06:02 UT, label "back" in Figure 1). We 
call the region between "outs" and "back" simply the back region of the MC. 
Finally, there is a region with weak, but coherent, field between the "back" 
and "end" boundaries. The physical origin of these regions is analyzed in 
Section 3.5. 

2.3. Orientation and extension of the cloud 

To facilitate the understanding of the MC properties, we define a system 
of coordinates linked to the cloud in which Zcioud is along the cloud axis 
(with B^, 

,cioud > 0). We define the latitude angle (^) between the ecliptic 
plane and the cloud axis, as well as the longitude angle {(p) between the 
projection of the axis on the ecliptic plane and the Earth-Sun direction 
(xgse) measured counterclockwise. Then, when ^=90° {9=-90°) the cloud 
axis is parallel (antiparallel) to zgse and it points to the ecliptic North 
(South). When ^=0° the cloud axis is on the ecliptic plane, ip=Q° being the 
case of the axial field pointing toward the Sun, and ip=90° ((^=270°) when 
it points to the terrestrial dusk (dawn). 

Since the speed of an MC is nearly in the Sun-Earth direction and it 
is much larger than the spacecraft speed (which can be supposed to be at 
rest during the cloud observing time), we assume a rectilinear spacecraft 
trajectory in the cloud frame. The trajectory defines a direction d (pointing 
toward the Sun); then, we define ycioud in the direction icioud x d and Xdoud 
completes the right-handed orthonormal base (xcioud; ydoudi -^cioud)- We also 
define the impact parameter, p, as the minimum distance from the spacecraft 
to the cloud axis. Then, we construct a rotation matrix from the GSE system 
to the cloud system and obtain the components of the observed magnetic 
field in the cloud coordinates: -Ba-^cioud, -By,cioud, -Bz,cioud- 

The local system of coordinates is especially useful when p is small com- 
pared to the MC radius (R). In particular, for p = and an MC described 
by a cylindrical magnetic configuration B{r) = Bz{r)z + B^{r)(p, we have 
^cioud = ^ and ycloud = when the spacecraft leaves the cloud. In this 
particular case, the magnetic field data will show: -Bj, cioud = 0, a large and 
coherent variation of -B^^cioud (with a change of sign), and an intermediate 
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and coherent variation of Bz^doud, from low values at one cloud edge, taking 
the largest value at its axis and returning to low values at the other edge. 

The minimum variance (MV) method (Sonnerup and Cahill, 1967) has 
been used to estimate the orientation of MCs (see e.g., Bothmer and Schwenn, 
1998, Lepping et al, 1990, Farrugia et al, 1999, Dasso et al, 2003, Gulisano 
et al, 2005). It gives a good estimation Hp is small compared to R and if 
the in-/out-bound magnetic fields are not significantly asymmetric. 

However, when the cloud presents a strong expansion, as in the event 
studied here, the directions derived by the MV method will mix two different 
effects in the variance of i?: (1) the effect of the coherent rotation of B (which 
provides the cloud orientation), and (2) the effect of the cloud 'aging' (the 
decrease of the field strength with time due to magnetic flux conservation 
combined with cloud expansion). This latter effect is not associated with the 
cloud orientation; thus, we apply the MV technique to B/B to decrease the 
cloud 'aging' consequences. 

We start the analysis taking the "in" and the "outs" boundaries, since 
the strongest magnetic discontinuity is located at "outs". With these bound- 
aries, we find the typical shape of the components of B in the cloud frame, as 
discussed above (see Section 3.3 for further justifications) . The MV method 
applied to the normalized field gives ^=-23° and (^=274°. With the same 
procedure and an end boundary on Nov. 10, lOrOOUT, Harra et al. (2006) 
found 6'=-20° and <^=276°. Changing the cloud end between "outi" and 
"outs", the ranges for 9 and (p are: 9 ~[-25,0]° and (p -[260,280]°. 

2.4. The data in the cloud frame 

Figure 2 shows the components of the magnetic and velocity field in the cloud 
frame for an orientation of the cloud axis such that: 9 = —10° and (p = 275° 
(this orientation is justified in Section 3.2). The magnetic field components 
show the typical large scale shape of MCs when the impact parameter is 
small compared with the cloud radius. However, an inner and non- typical 
sub-structure is present at its center (mainly observed in i?y,cio«<i) where it 
is antisymmetric). 

A strong expansion is observed in the velocity components, mainly along 
^^cioud- This is expected from the MC orientation and the data in GSE since, 
as if ^ 270°, the spacecraft cannot observe different parcels of fiuid along a 
large range of Zdoud values. However, a weak signature of expansion along 
the cloud axis can be observed in Vz,doud (Figure 2). For an MC with 9 ^ ^ 
and if > 270°, an axial expansion is characterized by Vz^cloud < in the 
front, changing to Vz^doud > in the back. This is just what is observed in 
our case when the mean value of Vz^doud is removed within the cloud. This 
gives us the clue that tp is larger than 270°. 
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Figure 2. Prom upper to lower panels, observed magnetic and velocity field components 
in the cloud frame. In this frame Scioud is along the cloud axis {0 = —10° and (p = 275° in 
GSE coordinates), ycioud is orthogonal to both the MC axis and the spacecraft trajectory, 
and icioud completes the right-handed orthogonal base. Vertical lines correspond to the 
same times as in Figure 1. The region between "in" and "center" is called the in-bound, 
and between "center" and "outi" or "out2" or "outs" the out-bound. Prom "outi" or 
"out2" or "outs" up to "back" we have the back region, while after that and up to "end" 
completes the ICME extension. 



When the spacecraft crosses a cyhndrical MC (or an elhptical one with 

one of the main axis parallel to the S un-Eartli direction) cuidp — 0, Bx^doud 
0. The first panel in Fig ure 2 shows that the observed Bx cloud 

has a slightly 

negative mean value. The sign of Bx^doud, together with the evolution of 
By^cloud, implies that the flux rope axis is above the ecliptic plane. 



3. Results with the Direct Method 

3.1. The Direct Method 

In this section we summarize and extend the direct method presented by 
Dasso et al. (2006). This method lets us find the rear boundary of a flux 
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rope for a given axis orientation, or the reverse, the MC orientation for 
a given position of the rear boundary. The front boundary of a flux rope 
is usually well defined by a discontinuity of the magnetic field (changing 
abruptly from a fluctuating field in the MC sheath to a strong and coherent 
field within the MC). The corresponding current sheet is expected to be 
present all around the flux rope as the limit between magnetic regions 
with different magnetic connectivities; so, with different magnetic stresses. A 
corresponding magnetic discontinuity, labeled "out" , is then expected at the 
rear of the MC. The flux rope is present in between these two discontinuities 
and the same amount of azimuthal flux is traversed twice by the spacecraft. 

Let us consider a flux rope at a given time. The conservation of the 
magnetic flux (V • B = 0) across a plane formed by the spacecraft trajectory 
and Zcioud {y cloud constant) gives (Dasso et al., 2006): 

/ -By, cloud dx dz = 0, (1) 

J flux rope 

with X, z being the spatial coordinates in the Xdoud and Zdoud directions, 
respectively. 

The observations provide only -B^^, cloud as a function of time along the 
trajectory. So, we need two hypothesis: an invariance of cloud along the 
flux rope axis and the conservation of the magnetic flux with time. The flrst 
hypothesis is justifled by a low ratio of the MC radius over the expected 
curvature radius of the axis and the balance of magnetic torques, which 
is expected to homogenize the fleld along the axial direction (Dasso et al., 
2006). The second hypothesis is valid as far as the amount of magnetic flux 
reconnected during the spacecraft crossing is low. Indeed, we have found that 
magnetic flux is reconnected in the front of the MC (sec Section 3.5). An 
estimation of the amount of flux reconnected during the spacecraft crossing 
gives ~ 5% of the initial azimuthal flux (see Section 3.6); so, unless the 
reconnection rate is much higher during the observing time than it was 
during the travel from the Sun, the amount of reconnected flux is small 
during the crossing. We neglect below such reconnected flux. 

The elementary flux crossed during dt is By^aoud{t) L{t)Vx^cloud{t) dt, 
where L(t) is the axial length of the portion of the flux rope which had a 
length Lin = L(tia) when the spacecraft entered the MC. Then, Equation (1) 
becomes: 

/ By^doud{t) L{t) Vx,cloud{t) dt = . (2) 

^flux rope 

If the axis orientation and the position of one boundary of an MC are 
known, the above flux balance property can be used to flnd the MC center 
and the other boundary as follows. We define the accumulative flux per unit 



dassoetal.tex; 1/02/2008; 11:48; p. 9 



10 



Dasso et al. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

Distance (AU) after:09-Nov-2004 20:30:19UT Distance (AU) after:09-Nov-2004 20:30:19UT 



Figure 3. Component of B perpendicular to both the trajectory of the spacecraft and to 
the cloud axis {By^cioud, thin line curve) and accumulated magnetic flux of this component 
per unit length (K^/Lin, thick line curve). The dotted vertical line marks the cloud center 
(Nov. 10, 02:02 UT). Left and right panels correspond to computation of F^/Lin without 
{L{t) = Lin) and with axial expansion (Equation 4), respectively. In both panels the cloud 
orientation was taken as ^ = —10° and tp = 275°. The vertical lines indicate the same 
positions as in Figures 1 and 2. 



length: 

~~? ~ / ^y,doud[t ) -J yx,doud(t ) Qt , [6) 

Lin J tin ^in 

where ti^ is the time of the MC front boundary (located at x = X^^) and x = 
It- Vx,cloudit') dt'. The position where Fy{x)/L\n has its absolute extreme 
gives an estimation of the position where the spacecraft reaches the closest 
distance to the MC axis. This indicates the x position of the MC center, 
being this estimation more precise as the impact parameter is lower. Then, 
when Fy [x) / Lin goes back to zero at a; = X^nt , we have the other boundary. 
The region from x = X[n to x = Xout defines the MC flux rope. 

Since the ratio of the MC radius to the Sun distance is small (typically ~ 
0.1 or lower), L{t')/L\n ~ 1- Below we derive a correction to this estimation. 
If the MC axis does not change drastically its shape, and if it does not 
disconnect from the Sun during the crossing time, its length is evolving 
proportionally to its distance to the Sun, D[t), so L{t)/Lin = D{t)/Din 
(where An ~ 1 AU for the present observations). Moreover, the cloud global 
velocity, which is the velocity of its center Vc, is not expected to change 
significantly during the crossing time; then, we have: 

All terms in the right hand side of Equation (3) can be derived from obser- 
vations. 
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3.2. Refined orientation of the cloud 

The direct method was previously applied to the Oct. 18-20, 1995, MC 
(Dasso et al, 2006). The orientation of its axis was well determined. The 
strong frontal discontinuity in the magnetic field was naturally related to an- 
other strong backward discontimiity, the flux balance given by Equation (2) 
was satisfied, and L(t) was constant. 

The orientation of the axis of the MC studied here is not so well deter- 
mined; then, we explore different orientations to find which angles {0 and ip) 
give a cancellation of Fy at the strongest magnetic discontinuity observed at 
the MC rear (labeled "outa" in Figures 1, 2), at the inner extreme labeled 
"outi"', and at the intermediate time "out2" (between "outi" and "outs"). 
Due to the orientation of the MC (i.e., its axis almost lying on the ecliptic and 
perpendicular to the Sun-Earth direction) Fy is mostly affected by the value 
of 6. We find 9 = —10° it 10° from variations of the end boundaries in the full 
range between "outi" and "outs", and using the two extreme possibilities 
on the axial expansion (no axial expansion, L(t) = Lin, and axial expansion 
proportional to the distance to the Sun, Equation 4). The value of ip is 
constrained by imposing that Bx^doudit) should have a small variation with 
time (as expected in flux rope models); so, with no contribution of azimuthal 
or axial field components. This gives ip = 275° it 10°. 

Left panel of Figure 3 shows Fy for L{t) = Lin (no axial expansion) 
and an orientation such that Fy is canceled at the discontinuity "outs", 
which gives 9 = —10° and ip = 275°. The rear boundary of the MC is 
in fact ambiguous from the data, since there is also a strong discontinuity 
between "outi" and "outa" ("out2", at Nov 10, at 09:00UT), see Figures 1 
and 2. Fixing the previous orientation, but using the axial expansion given 
in Equation (4) (see panel in Figure 3), we find that the cancellation of Fy 
is now at "out2" . With this boundary, we find a much less extended region 
with a reversal of cloud at the back of the MC. So, this boundary gives 
an MC field closer to the one inferred from a classical MC model. In this 
case, the results of the direct method agree with those of the fitted models 
(including expansion. Section 4). These are evidences that we have identified 
the right end boundary and orientation of the flux rope. 

3.3. Structure of the cloud 

The field component -Bj;,cioud vanishes at three locations near the cloud center 
(close to the vertical dotted line at the abscissa ~ 0.1 AU in Figure 3). Taking 
into account the expected antisymmetry of the in- and out-bound regions, 
we set the cloud center at Nov. 10 at 02:02 UT (dotted line in Figure 3). 
Then, the accumulated flux Fy/L\n (Equation 3) gives a unique relationship 
between the in- and out-bound data, since it labels each flux surface. This 
relationship is better shown using FyjL^a in the abscissa and reversing the 
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Figure 4- In-bound (red) and out-bound (blue) ydoud component of B in function of the 
accumulated magnetic flux of this component per unit length (Fy/Lin)- Fy/Li^ is computed 
as in Figure 3. The difference in the amplitude of By^cioud between the in-/out-bound 
regions is due to the MC aging (not corrected here). 

sign of By^cloud inside the out-bound branch (Figure 4). Then, peaks and 
valleys of i?y,cioud in the in-/out-bound branches can be related. As expected, 
this association is stronger near the MC center where the regions that are 
crossed are closer and also more isolated from the interaction with the solar 
wind environment (Dasso et al, 2005a). Except for the strong discontinuity 
at the MC borders, the association between structures is not clear outside the 
core (where cloud has no characteristic variations that can be recognized 
in both the in- and out-bound branches). 

There is a clearly distinguishable sub-structure in the cloud center seen 
as reversed peaks in -Bj^^doud and as a valley in Fy/L\^ (Figure 3). It has a 
small extension, ~ ±10"^ AU, and is globally anti-symmetric. A less evident 
sub-structure was also present in the previously analyzed MC (18 Oct. 1995, 
Dasso et al, 2006). 

li p = 0, the central sub-structure would imply the presence of a small 
twisted flux tube with opposite magnetic helicity in the center of the flux 
rope. The formation of such a structure is not possible in the corona. After 
analyzing several possibilities, the simplest interpretation is the following. 
Close to the minimum approach, the spacecraft trajectory is nearly tangent 
to the magnetic flux surfaces of the flux rope (-By^cioud ^ 0). Any warping 
of the flux surfaces gives a clear signal in the -By^cioud component. For geo- 
metrical reasons, such warping is more difficult to detect outside the center 
(where -By^cloud is important). It is noteworthy that such structure will be 
evident only if the magnetic data are rotated to the correct MC frame to 
have no mixing with the strong -B^, doud component. 

For the 18 Oct. MC, the warping was moderate and -By^doud kept its sign 
inside both the in- and the out-bound regions. For the Nov. 9-10 MC, the 
warping is more marked as sketched in Figure 5. This figure shows (to the 
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Distance (AU) after " in' boundary ■ 

Figure 5. Scheme of the magnetic field hnes near the MC central region (to the right) 
and the corresponding Fy/Lin evolution along the spacecraft trajectory (low left figure, 
a zoom of Figure 3). The magnetic structure of the fiux rope is undulated, probably due 
to its fast evolution and interaction with the surrounding medium. The undulations are 
amplified in this scheme for clarity, but in the observed MC they are significant enough 
to produce reversals of -Bj,, cloud- Dot at CC indicates the location of the cloud center. 
Numbers 1, 2, 3, 4 and 5 are reference points of the observed field lines. 

right) a scheme of the spacecraft trajectory (vertical dashed hne) across the 
MC core field lines (solid warped lines). The core crossing starts at field 
line '1' and ends at '5'. The large dot called 'CC marks the cloud center 
position. It is located towards positive values of Udoud since the cloud axis is 
above the ecliptic plane (see last paragraph of Section 2.4). Considering the 
cloud as a flux rope, which is compatible with the observations, the same 
accumulated flux Fy/Lin implies that the in- and out-bound field lines are 
connected (i.e., they are in fact the same field line observed twice). Thus, 
for example '1' is connected to '5'. The flux Fy/Lin has a local maximum 
at '2' and '4' and a local minimum at '3' when the spacecraft trajectory is 
tangent to the field lines. 

3.4. Magnetic fluxes 

The total azimuthal flux i<0 is estimated taking Fy{xcenter), assuming that 
p = 0. The largest source of uncertainty is the length of the flux rope, Lin, 
and in a more general way the assumed invariance by translation. Since this 
cloud presents signatures of being detached at one of its legs (Harra et ai, 
2006), we assume an initial length Lin = 1-5 AU, so intermediate between 
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values used in previous papers. The uncertainty in the MC boundary of the 
out-bound branch has a lower effect since the difference of found with 
the boundary "outi" and "outs" is only about 10% (see Table I). 

The axial flux F^, across a surface perpendicular to the cloud axis, can 
be estimated directly from the observations assuming a circular MC cross 
section, p = and neglecting the expansion to compute r = x(t) — Xcenter- 
Here, we neglect the axial flux in the core since it is a correction of the 
order of {p/Rf, see Dasso et al. (2006). Then, for the out-bound branch we 
compute F^^-^'^'"^ as (idem for but with the integration between 

tin ^nd ^center)' 



As in previously studied MCs (e.g. Mandrini et al. 2005, 2007; Attrill et al., 
2006) Fz is one order of magnitude lower than F^ (see Table I). 

3.5. Structure of the ICME 

The analysis of Fy/L\-a in Section 3.3 indicates that this MC is not formed 
only by a flux rope. Some of the MC characteristics (see Figure 1), such 
as a very unusually high magnetic field with a low variance, a low and 
a strong expansion, continue well after the rear boundary of the flux rope 
("outi", "out2", and "outs"). We discuss below the most plausible physical 
scenario to create such magnetic structure. 

There is an extended region where i?.y,cioud has still a negative and coher- 
ent behavior from position "outs" to "back", so the accumulated flux keeps 
increasing monotonously (Figure 3). This behavior was also found in the 
18 Oct. 1995 MC and it was interpreted as the trace of an original larger 
magnetic flux rope whose front was partially reconnected with the overtaken 
magnetic flux, as shown in Figure 6 of Dasso et al. (2006). In the example 
analyzed here, this interpretation has even more support from the data, 
since part of the overtaken flux is still present in front of the MC and we 
can estimate when reconnection started. 

There is a coherent negative -B^^ qIq^j field just in front of the MC (from 

18:00 UT to 20:30 UT, Figure 2). We interpret this as the remnant of the 
magnetic flux which reconnected with the original flux rope. At the time of 
observations, the MC overtakes this structure with a velocity difference of 

40 km/s; so reconnection is driven. Such velocity difference is crucial for 
reconnection efficiency since its rate increases with a larger velocity differ- 
ence (Schmidt and Cargill, 2003). Furthermore, Burlaga (1995) and Farrugia 
et al. (2001) have found that magnetic holes, such as the one preceding the 
MC analyzed here (top panel of Figure 1), are associated with magnetic 
reconnection. 
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Taking into account the previous reconnection scenario, the back region, 
where -By. cloud has still a coherent negative value, is simply the magnetic 
flux at the periphery of the original flux rope. This back region keeps the 
properties of typical MCs without fitting in the standard flux rope models. 
If the front reconnected, this back region is connected to the solar wind 
field; so, the magnetic field direction can change because of the propagation 
of Alfven waves, as observed here with the reversal of -Bz,cioud (Figure 2). 
However, if there is no other reconnection process, the flux of -B2y,cioud cannot 
be removed from the back of the MC; in fact, we find no other field to allow 
such reconnection. Thus, the closed (flux balanced) flux rope observed at 
1 AU is embedded in a larger structure which includes an extended back 
region. This flux rope was part of a larger one that was partially pealed at 
its front because it reconnected with its environment. 

3.6. Clues for magnetic reconnection 

The high velocity of the MC (Vc ~ 800 km/s), compared to its surroundings 
(Ri 600 km/s) implies a progressive extension of the back region (between 
boundaries "outs" and "back" in Figure 2). Assuming that the relative 
velocity, Vc — Vjaack) was similar at earlier times and that the -Bj^, cloud com- 
ponent (in the back region) had initially a value similar to that in the rear 
of the MC (because of pressure balance), we can estimate the period of 
time St between the start of reconnection and the MC observations. The 
range of time between the observation of the two extremes of the expanding 
region (between "outa" and "back") is Texpansion ~ 18 h. This corresponds 
to a spatial extension of ~ Texpansion Vback when the region was observed. 
So, the back region expands in size by Texpansion H)ack since the start of 
reconnection at the MC front. This extension has its origin in the relative 
velocity between the MC and the back region (which becomes lower because 
of its magnetic connection to the solar wind). From the beginning of the 
reconnection at the MC front, starting a lapse of time St earlier than the 
observations, the "back" boundary progressively separates from the MC with 
a relative velocity Vuc — Hack- Assuming that this relative velocity was not 
changing drastically during the MC transit, the back region expanded by 
~ (^MC — ^back)^i- Equating the previous estimations of the back region 
extension, we find: 

St ^ '"^"P""'"';^^"^'^ ^ 18 x 600/200 !^ 54 h . (6) 

Vc — Vback 

The solar event that is the most probable source for the MC studied here 
occurred in AR 10696. It is a multiple event that starts with a steep rise 
in GOES light curve at ~ 15:50 UT on Nov. 7, and has two clear peaks, 
one at ~ 16:00 UT and the other at ~ 16:35 UT. At the time of the second 
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peak a large two-ribbon flare was observed within AR 10696, but also two 
Ha ribbons were seen at both sides of an erupting trans-equatorial filament 
extending from AR 10696 to AR 10695 at the southwest (see Harra et al., 
2006). During this intense event that reached class X2.0 in soft X-rays the 
full neutral line, which formed a switch-back, erupted. The CME observed 
in the Large Angle and Spectroscopic Coronagraph (SOHO/LASCO) C2 on 
the Nov. 7 at 16:54 UT is associated with this multiple flare involving AR 
10696 and the erupting trans-equatorial filament to the southwest. Longcope 
et al. (2007) have proposed that the source of the MC is the AR eruption, 
while Harra et al. (2006) have considered both possibilities either the AR or 
the trans-equatorial filament eruption. We discuss below these scenarios in 
view of our IP analysis. 

Taking ~ 15:50 UT on Nov. 7, as the start time of the solar event, and 
an arrival time for the MC front at ~ 20:30 UT on Nov. 9, wc obtain a 
transit time of ~ 52 h, comparable to the transit time 47 h) computed 
assuming a constant MC velocity of -800 km/s and the spacecraft located at 
LI. Both estimations are comparable to 5t (Equation 6); this implies that 
reconnection started when the fiux rope was close to the Sun. 

What could be the magnetic structure present in front of the flux rope, 
probably already from its origin in the corona? The field in front of the MC 
is oriented mostly southward with an average field intensity of ^ -30 nT, 
though there is a structure of northward oriented field with a temporal length 
between 19:54 UT and 20:01 UT on Nov. 9, an average field strength of ~ 
-15 nT, and a spatial extension of 0.002 AU (taking an average velocity 
of ~ -800 km/s). Analyzing the MDI magnetogram on Nov. 7, the most 
probable origin of the southward oriented structure is the large scale, nearly 
potential field of AR 10696. If the source of the MC is the AR eruption; 
then, the core of the AR should become kink unstable (as it was observed 
in an eruption which occurred in the same AR three days later (Williams 
et al, 2005)). In the kink instability, part of the twist is transformed into 
writhe implying a strong rotation of the flux rope. Taking into account the 
MC orientation found in Section 3.2, this rotation should be ~ 160°. The 
MHD simulation of Gibson et al. (2004) gives a writhing of the flux tube of ~ 
120° (before it reconnects with the overlying fleld). Numerical simulations 
by Torok and Kliem (2005) confirm this using a different approach and, 
indeed, the writhing could be as large as a rotation of 160° depending on 
the properties of the overlying field (Torok, private communication) . So the 
amount of rotation and its direction (for a left-handed fiux rope), if the MC 
fiux rope comes from the AR, is coherent with recent MHD simulations of 
kink unstable flux ropes. 

Interplanetary scintillation observations suggest that the material from 
the core of the AR was ejected primarily northward and, thus, it could 
remain unobserved at 1 AU (see Harra et al., 2006). In this scenario, and if 
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the source of the MC is the trans-equatorial filament eruption, the direction 
of the trans-equatorial filament and the cloud axis differ by 67*^ (see Figure 
12 in Harra et al., 2006). Moreover, the sense of rotation from the filament 
to the cloud axis is opposite to the one expected from the development of a 
kink instability in a region of left-handed magnetic helicity. As discussed by 
Harra et al. (2006), it may happen that the kink instability has not played 
a role far away from the strong fields of the AR. The northward oriented 
structure observed between 19:54 UT and 20:01 UT on Nov. 9 could be 
related to the northward oriented trans-equatorial loops observed above the 
filament. However, in both scenarios rcconncction should be forced between 
the ejected flux rope and the above AR arcade field from the beginning of the 
eruption. This explains the origin of a nearly antiparallel field (neglecting 
the short period of northward directed field) in front of the flux rope from 
the beginning of the launch, as indicated by a value of 6t similar to the 
transit time. 

Our previous discussion has implications for the reconnection rate in a 
collisionless plasma. While the erupting flux rope pushed against an overly- 
ing nearly anti-parallel magnetic field, we still observe part of this overlying 
flux in front of the MC at 1 AU! Of course we have no way to determine 
how strong was the forcing, nor the time dependence of the reconnection 
rate during the flux rope transit. Still this observation is a clue that mag- 
netic reconnection is not efficient in a collisionless plasma, as expected from 
classical theory. Some recent observations show direct evidence of magnetic 
reconnection in a collisionless plasma as the solar wind (Gosling et al, 
2005). Some numerical simulations show that the Hall effect can increase 
the reconnection rate above the classical rate (e.g.. Morales et al., 2005). 
Our observations set a constrain on the reconnection rate, which can be 
quantified only after a numerical modeling of the flux rope ejection and 
transit to 1 AU. 

Prom the computed accumulated flux (Figure 3) we can estimate the 

relative amount of reconnected flux from the original flux rope ejected from 
the Sun as « 1.25/7.4 17%, assuming L{t) = (left panel). With 
the expansion length given by Equation (4), the flux present in the back 
part, which is the fraction of reconnected flux, is ~ 30% (however it is not 
obvious that we can use Equation (4) for the back part where the field is 
connected to the solar wind) . The relative amount of flux present in front of 
the MC is ^ 0.6/7.4 8%. This implies that the relative amount of stable 
(not kinked) magnetic flux, arch-like, above the erupting flux rope is in the 
interval [25,38]% of the azimuthal flux in the MC. 

Finally, a fraction of the -Bj^,cioud magnetic flux observed after the bound- 
ary "in" is expected to be reconnected when the spacecraft exits the flux 
rope. To estimate an order of magnitude for this flux, let us suppose that 
the reconnection rate was comparable to its mean value during the transit 
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from the Sun. The amount of reconnected flux during the MC observing 
time is the ratio of the crossing time (~ 14 h) over the transit time from the 
Sun (~ 52 h) times the flux reconnected from the Sun. Then, ~ 5% of the 
original flux was reconnected during the MC crossing time. So, the fraction 
of flux reconnected between t = tin and tout is expected to be an amount 
much lower than the uncertainties in the flux estimations (see Table I). 



This cloud presents a large velocity difference between its front and rear 
parts (fourth panel in Figure 1). This velocity difference is the consequence 
of a large MC size and a significant expansion. Both conditions imply that 
the magnetic field is observed at significantly different times during the 
MC crossing, so at times when the MC has significant different sizes. As a 
consequence of magnetic flux conservation, this implies the observed decay of 
B with time (upper panel of Figure 1). Then, Bz^doudi and even more By^^loud 
(since it is stronger close to the flux rope borders) presents a remarkable 
asymmetry between the in- and out-bound branches. Moreover, the center of 
the cloud is observed before the central observing time for the full structure, 
as expected for an spatially symmetric expanding object. 

In this section we compare the observations to fitted models describing the 
evolution of the magnetic field assuming an isotropic self-similar expansion. 
The observed velocity is Tiscd to derive the expansion rate, which is then 
used in the expansion magnetic model. We also quantify the magnetic fluxes 
from the fitted models, and compare them with the results obtained from a 
classical static model and the direct method. 

4.1. Expansion Model 

We assume an isotropic self-similar expansion of the MC, where all distances 
are multiplied by a factor f{t): 



where /(to) = 1 and ro is the position of a given element of fluid at some 
reference time tg- Since the MC front boundary, called "in", is well defined, 
we select this as the reference time, to = ^in- Each element of fluid is labeled 
by its rb value, then the ro coordinate is a marker of each element of fluid 
(Lagrangian coordinates). 

We assume below that the flux rope size increases linearly with time. 
Then, /(t) can be written as: 



4. Results using fitted models 



r(t) = ro/(t), 



(7) 



/(t) = 1 + (t - tin)/r 



(8) 



dassoetal.tex; 1/02/2008; 11:48; p. 18 



Transformation of a flux rope to an ICME 



19 



where T is the time of expansion. If Equation (8) would be valid for all 
the transit time from the Sun, then the flux rope size would vanish at time 
i = ^in — T, so T would also be approximately the transit time from the Sun 
to the position observed at t\^. In practice, we use Equation (8) only during 
the MC crossing; so, we are only assuming a local linear increase of the size 
with time. Depending on the evolution during the transit from the Sun, T 
could be diflFerent from the transit time. 

The time evolution for the plasma velocity V in the cloud frame comes 
from Equation (7) doing the temporal derivative and keeping ro fixed (so 
following a plasma element): 

^(0 = 7^ — = ^ • (9) 

The first equality gives the velocity at a given time t, so it is proportional to 
the distance. The velocity decays with time, but this decay is only apparent 
since when one follows a given plasma element, defined by fg, the velocity 
is in fact constant (second equality). Then, Equation (8) implies that there 
are no forces acting on any plasma element (free expansion). This is an "a 
posteriori" justification of Equation (7); for an isotropic expansion, an initial 
force-free field stays force-free. 

The observed speed is the sum of the expansion speed and of the global 
speed of the MC. During the time of observation, we can assume that the 
MC is globally moving at a constant speed — K (velocity of the MC center) 
along Xcioud because the aerodynamic drag significantly affects Vc only in 
a time scale comparable or larger than the transit time to 1 AU. We also 
assume that the spacecraft is at rest. Then, the velocity component along 

•^cloud 

This equation is fitted to the data to derive Vc, tc (the time when the 
spacecraft crosses the cloud center) , and T. It is noteworthy that the nearly 
linear observed Vx is a consequence of the constant speed Vc, and not of the 
linear expansion with time assumed in Equation (8). The assumed expansion 
introduces only a non linear correction in (t — tin) because the crossing time, 
tout — tin, of the MC is small compared to the time of expansion, T (a further 
correction would come from the non-linear development of a general f{t)). 
This justifies the use of Equation (8) 

4.2. Results for the expansion 

The observed Vx^cioud has a globally linear variation with time within the MC 
(Figure 6). To better show its variation, we subtract the mean value of the ve- 
locity computed inside the range of positions "in" to "outs", 
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Figure 6. Observed (dots) and fitted (thick dasli-dotted line) radial velocity profile 
{Vx,cioud— < Vx,cioud >)• Vertical lines mark the same times as in Figure 1. 



< y-xjdoud >= —794 km/s. When the impact parameter (p/R) is small, 
as it is in this MC, this speed represents the radial velocity with respect to 
the cloud axis, such that for a cloud in expansion, Vx^cioud— < yx,cioud > is 
negative before the spacecraft reaches the cloud axis and positive after that. 

From a least square fitting of Equation (10) to the velocity data inside 
the "in" to "outs" boundaries, we obtain T « = 73 h (rs 3 days). We find 
a slightly longer time of expansion T Ti = 79 h (pa 3.3 days) when the fit 
is restricted to the interval between the "in" and "outi" boundaries, and we 
get T2 = 77 h (~ 3.2 days) for boundaries "in" to "out2". Considering the 
self-similar expansion given by Equation (8), this implies that the MC has 
expanded from the time tin by a factor 1.19 and 1.15 when the spacecraft 
crossed "outs" ^^'^ "outi", respectively. 

Both times of expansion, and Ti, are longer than the transit time from 
the Sun (si 52 h. Section 3.6). Taking an approximately constant global 
velocity Vc, the distance of the MC front to the Sun increases by the factor 
D(t) /I?in ~ 1 + (t — ti^)Vc/ (see Equation 4). Doing the ratio of the mea- 
sured expansion rate variation (^^^) to '^^f^ , we obtain an undimensional 
factor Din/{VcT). This factor is 0.64 and 0.59 for boundaries "outs" and 
"outi", respectively. Considering the orientation of the MC {6 = —10° and 
If = 275°), the velocity in Figure 6 is measured mainly across the fiux rope. 
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Then, the MC is expanding at a significantly smaller rate radially than what 
we expect assuming the isotropic expansion in Equation 4. 

Assuming that the spacecraft trajectory is close enough to the MC axis, 
p/R « 1, the observed velocity lets us compute the radius of the flux rope 
since the distance crossed (/ V;j;,cioud dt) is i?in + -Rout = -Rin(l + /(^out))- 
Taking the boundary "outs", we find i?in = 0.12 AU and i?out3 = 0.14 AU; 
while for boundary "outi" we find smaller radii, R\n = 0.10 AU and -Routi = 
0.12 AU. If we do not consider the expansion, we get the mean values because 
the distance crossed is simply 2R. 

4.3. Magnetic field models 

Lundquist's static model (Lundquist, 1950) is a classical linear force free 
configuration (V x 5 = aB, with a constant). From the conservation of 
the magnetic flux during the expansion, and assuming a Lundquist's field 
at a given time, it is possible to derive an expansion field model (see, e.g., 
Shimazu and Vandas, 2002, Berdichevsky et al, 2003). We will also consider 
a more general model where the amplitude of the azimuthal and axial com- 
ponents are independent. This represents an approximation for a flux rope 
with an oblate cross section (see Vandas and Romashets, 2003 for an exact 
solution). We consider this modified expansion Lundquist's model to keep 
the same functional dependence of the field components, so the difference 
between the three models is only the number of free parameters (and the 
physics involved). 

The three models are described by the equations: 



where -Bin,(/,, 5in,z; and are the field and a values when the spacecraft 
enters the MC at t = ti^. For the static model f = 1 and i?in.</, = -Bin.z = Sin, 
so there are only two free parameters, Sin and ain. For both expansion 
models f = f{t) is function of time as given by Equation (8). For the 
expansion Lundquist's model, -Bin,,}!) = -Bin,z = -Bin, so there are also two 
free parameters as for the static model. For the modified model, there is an 
extra free parameter; the parameters are: -Bin,,/,, -Bin,z and ain- Notice that 
for all models -B^ is not forced to vanish at the MC boundaries (this lets ccin 
as a free parameter). 

Because for a fixed time t, each of the components and B^ have a 
spatial dependence as in Lundquist's model, the equations for the magnetic 
fluxes Fz and F^f^ are the same as in Dasso et al. (2006), but now R = R{t) = 
Rin f{t) and L = L{t) = Lin f{t) with R^ and -Lin the radius and the length 



Br{r,t) 
B,{r,t) 



0, 

-Bin,,/, Ji{ahi r/i) 
Bin,z ^o(ain r/i) 



(11) 
(12) 

(13) 
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Figure 7. By^cioud (left panels) and B^^cioud (right panels). Observations (dots), cylindrical 
dynamical model (red dotted lines), cylindrical static model (blue dotted lines), and 
modified model (see text, green dotted lines). Upper panels show the fitting using the 
end time as Nov. 10, 10:20 UT (outs), lower panels use Nov. 10, 8:15 UT (outi) as the 
end time. Vertical lines mark the same times as in Figure 1. 



of the cylinder at time t = tin- Then, Fz and are, as expected, constants 
of motion, because the increase in L{t) and R{t) cancels the decay of the 
field components and of a. They both simply write: 



27rBin,0i?in^l(ain-Rin) 

= , (-14) 

Bin z (1 — Jo(ain.Rin)) ^ /i r-\ 

F<p = Lm. (15) 



4.4. Results for the magnetic field 



In this section we fit the MC observations with the three models described in 
Section 4.3. We fix the orientation to the one given in Section 3.2 (0 = —10° 
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and if = 275°). This allows us to test the effect of including the expansion 
and the decoupling of azimuthal/axial field from the problem of finding the 
MC axis. 

We first use the data in between "in" and "outs". We use a nonlinear 
fitting routine to fit the models presented in Section 4.3, assuming p = 0, to 
the observations of By^cioud (which corresponds to ±B^^cioud) and Bz^doud- 
The static model cannot reproduce the observed asymmetry due to the decay 
of the field and the shift of the position of the cloud center, because of its 
intrinsic symmetry (top panels of Figure 7) . Both the static and expansion 
Lundquist's models overestimate the axial field, Bz.doud: near the cloud 
center and underestimate the azimuthal field, By^cloudj near the boundaries. 
When the extra freedom -Bin,z 7^ ^in.<^ is considered, the model can reproduce 
significantly better the observations. 

A quantitative comparison between the different models is given by \/^, 
where is the time average of (B^odel-Bobservations)'^ ■ We find V3?=14:16:19 nT, 
for the modified model, expansion and static Lundquist's models, respec- 
tively. For the three models a shift of the positions where By^cloud = and 
where B^^doud is maximum to a time later than in the observations is present 
(top panels of Figure 7) . 

We also explore the two earlier rear boundaries: "outi" and "out2". A 
significant better fit is found with the boundary at positions "outi" and 
"out2" (lower panels in Figure 7 show the fitting for "outi", very similar 
fitted curves are obtained for "out2"). 

A quantitative comparison between the fits is given in the form 

^ Xout,i • \/ Xout,2 • \J Xout,3 i^ units of nT. For the static model we obtain 
16:17:19, for the expansion model 15:15:16, and for the modified model 
12:12:14. 

If the boundary is located at "out2", the modeled cloud center corre- 
sponds now to Nov. 10 at 02:16UT, 14 minutes later than the center given 
by the observations and the direct method in Section 3.3 (dotted line, Nov 
10, 02:02UT). For a rear boundary at "outi", the modeled cloud center 
corresponds to Nov. 10 at 01:58UT, only 4 minutes earlier than the value 
obtained using the direct method. 

Let us now analyze the differences between the observations and the best 
model (modified Lundquist's model with boundary "outi"). The model fol- 
lows globally well the observations, except for the central reversal of By^doud 
and close to the boundaries. This central reversal cannot be taken into ac- 
count by the model (see Section 3.3 for an analysis of this feature). Close to 
the boundaries, the asymmetry of the model is not as large as the observed 
one; the isotropic expansion model gives By^oMti/ By^i^ = /"^(touti) ~ 0.75, 
while the observations give « 0.5. Due to the crossing geometry, the data 
mainly reflect the radial expansion velocity. Let us now consider a refinement 
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Table I. Magnetic fluxes present in the flux rope in units of 10 Mx using 

different fitted models or the direct method, for three backwards boundaries 
("outi", "out2", and "outs" in Figures 1 and 2), and for orientation given 
hy 6 — —10° and ip — 275°. For the direct method we present an average 
between the values obtained for the in-bound and out-bound branches {F^ is 
not affected by including axial expansion). The fluxes are not corrected for the 
flux lost by magnetic reconnection during the MC travel from the Sun (since 
only the direct method permits an estimation of the reconnected flux) . The 
three models are described by the same equations but differ by the constraint 
set on the free parameters, so the physics involved (Equations 11-13). An 
initial length for the cloud Lin = 1.5 AU is assumed. 



end: "outi" "out2" "outs" 



Model or Method 






F, 




F, 


F^ 


static 


7.4 


60 


7.7 


60 


7.4 


60 


expanding 


7.4 


64 


7.4 


67 


6.8 


69 


modified 


6.4 


91 


6.8 


96 


6.2 


100 


direct without axial expansion 


5.0 


81 


4.7 


85 


2.1 


91 


direct with axial expansion 


5.0 


90 


4.7 


95 


2.1 


102 



of the above model, the radial expansion is still given by f{t), while the 
axial expansion is rather given by g{t) = L{t)/L^^ (Equation (4)). By^cioud is 
affected both by the radial and the axial expansion. The conservation of the 
azimuthal flux gives: By{t) = By^\ii/[f{t)g{t)]. In the fitting of the models, 
we have assumed f{t) = g{t), but in fact D\a/Vc < T so that g is slightly 
larger than /. With An/K ~ 47 h and T 79 h, we find /outi = 1-15 
and (7outi = 1-25, which gives i?y,outi/-Bj,,in ^ 0.69 (rather than 0.75 with an 
isotropic expansion), a value closer to the observed value (« 0.5) but still 
larger. 

A spatial asymmetry between the in- and out-bound branches, which 
cannot be attributed to the expansion, is present in the observations (Fig- 
ure 7). Indeed, the observed By^c.loud does not have the expected decrease 
towards the front boundary (which is present in the model and also in the 
observations towards "outi", Figure 7). At the MC front the magnetic field 
is expected to be compressed (then, it is enhanced) by the dynamic pressure 
of the overtaken plasma. The interaction with the surroundings is likely to 
be at the origin of this extra asymmetry (on top of the expansion) between 
the in- and out-bound branches. 

From the fitted parameters and expressions given in Equations (14)- (15), 
we obtain the values for the magnetic fluxes using the different models 
(Table I). We assume an initial length Ljn = 1.5 AU as in Section 3.4. 
From Table I, is estimated in the range [2.1-7.7] xlO^'' Mx, and in 
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the range [60-102] xlO^o Mx. Fixing the end boundary at Nov 10, 10:00 UT 
(only 20 minutes earlier than boundary "outs"), and using the static classical 
Lundquist's model, Longcope et al. (2007) reported = 7.2 x 10^° Mx and 
F^/L = 41 X 10^0 Mx/AU (thus F^ = 62 x 10^° Mx, for a cloud length of 
L = 1.5 AU, as assumed here). These values, as expected, are very close to 
our present results with the static model and boundary "outs" • Our present 
results show that the expansion affects slightly the computed fluxes (more 
F(p), while decoupling the fits of -Bj/,cioud and -B^,cioud has the largest effect. 
For rear boundaries "outi" and "out2", the estimations of both fluxes using 
the modified Lundquist's model arc in close agreement with the results of the 
direct method that consider an axial expansion, in particular this agreement 
is much better for F^. Considering the rear boundary "outs", the value of 
Fz obtained from the direct method is lower than the one from the modified 
Lundquist's model due to the significantly negative value of cloud beyond 
"out2" (see Figure 2); this contributes to decrease F^ when integrating the 
circular cloud section between "out2" to "outs". As discussed in Section 3.2, 
we believe that the rear boundary of the cloud should be between "outi" 
and "out2". 



5. Summary and Conclusions 

The ICME of Nov 9-10, 2004, was a complex event with a large expansion 
and a strong magnetic field in the front decreasing monotonously (almost 
linearly) with time. Earlier analysis considered that two magnetic clouds 
(MCs) were located inside this ICME. Latter studies concluded that only 
one MC was present. We confirm this and precise the orientation and bound- 
aries of the flux rope using several methods. We also find clues about the 
interaction of this MC with its surroundings. 

To facilitate the understanding of the physics involved in the MC it 
is useful to transform the data to the local MC frame where the axial 
and azimuthal components of the magnetic field are decoupled. A classical 
method to determine the MC axis is the minimum variance (MV) method, 
which takes into account the different spatial behavior of the magnetic field 
components to find the flux rope orientation. We minimize the effect of the 
strong expansion, which implies a decreasing field magnitude with time, 
normalizing the field at each data point. This gives a range of possible 
orientations (typically with a precision of the order of ±20°). 

The determination of the orientation of the flux rope axis is improved 
using the direct method. This method is based on two main points: first, 
the flux rope is topologically distinct from the surroundings, so it should 
generically be bounded by a discontinuity of the magnetic fleld components 
(presence of a current sheet), and second, the same amount of azimuthal 
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magnetic flux should be present in the in- and out-bound branches of the 
cloud. For the cloud of Nov. 9-10, 2004, the frontal discontinuity is well 
defined, while three rear discontinuities are present, called "outi", "out2", 
and "outs". They are separated by about 2 h compared to a MC duration 
of about 14 h. The azimuthal flux relates the frontal discontinuity to two 
discontinuities at the rear using two different hypotheses for the axial evolu- 
tion: an expansion comparable to the radial one gives a rear discontinuity at 
"out2" , and a negligible axial expansion gives a rear discontinuity at "outs" . 
However, since the second case implies a reversal of the axial held at the rear 
of the MC, we conclude that the first discontinuity ("out2") is associated to 
the frontal discontinuity. This defines precisely the extension, as well as the 
orientation angle 6 (within ±5°, ^ = 10 it 5°) of the flux rope. 

Fitting a model to the data is another approach to understand the ob- 
served magnetic structure. Wc have used three models that arc based on the 
Lundquist's solution. The first one is the classical static solution. The second 
one includes a self-similar expansion with the same rate in the axial and 
radial directions. Finally, the third one also includes an isotropic expansion 
and decouples the fit of the azimuthal and axial field components to take into 
account the observed stronger azimuthal component (a possible signature of 
a fiat cross section). The expansion rate is obtained fitting the model to 
the observed plasma velocity. The best fits to the data arc obtained when 
the first and second discontinuities (labeled "outi" and "out2") are used, in 
agreement with the results obtained with the direct method. 

Comparing the results of the fitted models with the direct method, using 
boundaries "outi" and "out2", we find that the axial and azimuthal fluxes 
are in the ranges [4.7-7.7] xlO^o Mx and [60-95] xlO^o Mx, respectively. The 
main limitation on the axial flux measurements is the unknown shape of the 
cross section. For the azimuthal flux, it is important to consider the axial 
expansion. Here the limitations are different, the shape of the cross section 
is not important, the main limitation is the distribution of the flux along 
the MC axis. Finally, we confirm that the azimuthal flux is one order of 
magnitude larger than the axial flux. 

After the large and coherent rotation of B, some typical MC character- 
istics are still present: low level of fluctuations, strong expansion (observed 
in the decay of V and -B), intensity of magnetic field higher than the typical 
solar wind values, low /3p. These are evidences that the MC extends farther 
in the back of the flux rope. Part of this back (from "outs" to "back" in 
Fi gures 1 and 2) shows a coherent behavior of By (^iqu^, which we interpret as 
the signature of an originally larger flux rope that was partially reconnected 
in its front near the Sun, with the consequent flux removal. 

Reconnection of the cloud field with the overtaken solar wind held is 
another source of underestimation of the original magnetic flux that was 
launched from the Sun. We have found that ~ 17 % (with a 30 % as an 
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upper bound) of the azimuthal flux was lost in the front of the MC during its 
travel from the Sun; this is much less than in the Oct. 18-30, 1995 MC, where 
it was estimated to be about 57 % (Dasso et al., 2006). Another difference 
between the MCs is that the Nov. 9-10 MC has a back part moving at a 
speed significantly lower (by f« 200 km/s) than the flux rope (while the 
Oct. 18-20 MC was overtaken by a fast stream). Then, the reconnected held 
progressively forms an extended region in the back of the flux rope (with 
a weaker magnitude but still with a smooth spatial variation). Prom the 
extension and velocity difference of this region with the flux rope, we have 
estimated that reconnection started close to the Sun, possibly between the 
erupting twisted flux tube (giving the flux rope) and the overlying arcade of 
active region 10696. Part of the arcade fleld is probably present in front of the 
MC, with a nearly anti-parallcl direction and a signiflcant velocity difference 
(~ 40 km/s), indicating that magnetic reconnection is not so efficient in the 
interplanetary space. 

Reconnection in front of the MC has several observational consequences, 
as follows: 

First, it introduces an asymmetry in the observed magnetic fleld. The 
remaining part of the flux rope is observed first followed by an extended 
tail of weaker magnetic field (which is re-orientated since it has changed 
its connectivity). When this process is dominantly at work, this implies the 
presence of a closed flux rope at the beginning of the ICME. 

Second, depending on the solar launch direction, the spacecraft could 
cross the flux rope or its laterally extended back part. In this last case, one 
would detect some characteristics of the MC (such as a coherent field and low 
/? values), but without the coherent rotation of the field. Such crossing would 
be classified as an ICME (without MC). An example of such observations 
with the two HELIOS spacecrafts is analyzed by Cane et al. (1997) and 
another example using the ACE spacecraft is analyzed by FouUon et al. 
(2007). 

Finally, the back flux connected to the solar wind field makes the mov- 
ing magnetic structure larger in the transverse direction (orthogonal to the 
global motion). Prom its mixed origin the back region is expected to move at 
a speed intermediate between the MC and the solar wind speed, as observed 
in the present analyzed MC. Then, with a significant larger velocity than 
the surrounding medium, the back region is expected to have an effect on 
the frontal shock surrounding the ICME. A larger transverse scale implies 
a forward shock at a larger distance in front of the MC than the distance 
deduced by its flux rope transverse size (Farris and Russell, 1994). Such 
large distance has so far being interpreted as a flat flux rope (Russell and 
Mulligan, 2002). The consequence of reconnection in the front of the MC is 
an alternative and / or complementary explanation which requires numerical 
simulations to be quantifled. 
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The results obtained show the potentiahty of combining several methods 
of analysis, minimum variance, direct method and fit to the data. This 
analysis will be done for other MCs to derive the variety of possible physical 
scenarios and also to improve our understanding of MCs and ICMEs. 
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